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Abstract 

Kinetic equations are often appropriate to model the energy den- 
sity of high frequency waves propagating in highly heterogeneous me- 
dia. The limitations of the kinetic model are quantified by the sta- 
tistical instability of the wave energy density, i.e., by its sensitivity 
to changes in the realization of the underlying heterogeneous medium 
modeled as a random medium. In the simplified Ito-Schrodinger regime 
of wave propagation, we obtain optimal estimates for the statistical 
instability of the wave energy density for different configurations of the 
source terms and the domains over which the energy density is mea- 
sured. We show that the energy density is asymptotically statistically 
stable (self-averaging) in many configurations. In the case of highly 
localized source terms, we obtain an explicit asymptotic expression for 
the scintillation function in the high frequency limit. 
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1 Introduction 



Let us consider the following scalar wave equation for the pressure potential 
p(r,x,t): 



where r is time, (x, t) eR^xR denote the spatial variables, A x is the Laplace 
operator in the transverse variables x, and c(x, t) is the local sound speed. 
Our objective is to understand the properties of p(r, x, t) when c(x, t) is a 
highly oscillatory random field and the initial conditions for p(r, x, t) oscillate 
at the same frequency. 

The analysis of high frequency waves in random media based on (JTJ) is 
extremely complicated and still not totally established mathematically. Since 
the wave field is oscillatory, its (weak) limit typically misses most of the 
energy of the wave field p. Kinetic models are then used to capture the 
energy density of the wave fields; see e.g. [HI El [12], EJ for rigorous results, 
[3] [19] for more formal derivations, and [HI [TU [THl [20] for references in the 
physical literature. 

The validity of the kinetic model is limited by its statistical instabil- 
ity, namely by its variability when the realization of the underlying random 
medium is changed. In many situations, the energy density is self-averaging 
[21 El El, which means that the energy density measured (averaged) on a 
(sufficiently large) domain is asymptotically, as the frequency goes to infin- 
ity, independent of the realization of the random medium. The above results 
often require that the domain of measurements be of size independent of the 
wavelength and that the source term for the kinetic model be sufficiently 
smooth. 

In this paper, we are interested in the statistical stability of such kinetic 
models in a very simplified regime of wave propagation, namely the Ito- 
Schrodinger regime. The latter regime arises when the wave field is a very 
narrow beam propagating in the direction t and the sound speed c(x, t) os- 
cillates more rapidly in the direction t than it does in other directions. Such 
assumptions are valid in somewhat restrictive practical settings. However, 
this regime of wave propagation is relatively simple to analyze mathemat- 
ically and provides interesting qualitative answers regarding the statistical 
stability of more general kinetic models. 

The validity of kinetic models has been analyzed numerically in several 
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settings [HI El ED], with quite good agreements with the energy density given 
by wave equations of the form ([T]). Such kinetic models may then be used to 
solve inverse problems, where constitutive parameters in the transport equa- 
tion modeling e.g. buried inclusions or statistics of the random medium, are 
reconstructed from available boundary measurements. We refer the reader 
to 0, [TU] for reconstructions based on synthetic (numerical) data and to [3] 
for kinetic reconstructions from experimental data in the micro-wave regime; 
see also [I] for a review on the use of kinetic models in the imaging of buried 
inclusions. These studies show that the kinetic models perform relatively 
well. Their limitations are almost entirely caused by our lack of knowledge 
of the random medium, which generates some statistical instabilities in the 
measurements. Understanding these instabilities will allow us to improve 
on the reconstructions and to have a better understanding of the maximal 
resolution that can be achieved. 

Ito-Schrodinger regime. In the Ito-Schrodinger regime, we introduce 
^(x, t; k) as 

P (t, x, t) = ±- / e^ c ^V(x, *; «)cbd«, (2) 

where Co is the background sound speed, assumed to be constant. Thus if) 
represents waves at position (x, t) propagating with frequency u> = cq\k\. 
After appropriate scalings and simplifications, the wave field tj) satisfies the 
following Ito-Schrodinger stochastic partial differential equation: 

d^uCM;*;) = - («77 A x - K 2 R(0))ip v dt + iKip v B(-,dt\. (3) 

Since k plays no significant role in the sequel, we set it to k = \. Here, 
£?(x, dt) is the standard Wiener measure, whose statistics are described by 

E{£?(x, t)B(y, t')} = R(x - y)t A t', (4) 

where E is mathematical expectation with respect to the measure of an ab- 
stract probability space on which S(x, dt) is defined and t A t' = min(£, t'). 
We shall not justify ([3]) from ([1]). See [1] for a justification in one dimension 
of space and [2] for the scaling arguments leading to 

For our purposes, i/j v {x.,t) satisfies a wave equation with highly oscilla- 
tory coefficients oscillating at a frequency inversely proportional to the small 
parameter r\ <C 1. We assume that ^(x, 0) also oscillates at a frequency 
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comparable to r]~ l and are interested in the properties of the wave field as 
rj — > 0. Because the field oscillates rapidly, its weak limit is of little interest. 
A more interesting quantity is the energy density of the waves \x/j v \ 2 (x,t), 
or the probability density in the context of quantum waves. Because the 
energy density does not satisfy a closed-form equation, it is more convenient 
to analyze energy densities by introducing the following Wigner transform of 
the wave field: 

W v (t, x, k) = JL J e**^ (x - f)^(x + "f, t) dy, (5) 

where i/j v denotes complex conjugation of ijj. Note that f Rd W^t, x, k)c£k = 
^(x, t)\ 2 by inverse Fourier transform so that may be seen as a phase 
space (microlocal) decomposition of the energy density. 

Let i/),(x, 0) be a sequence of functions uniformly bounded in L 2 (R d ), 
77-oscillatory, and compact at infinity in the sense of [13], i.e., such that for 
every continuous compactly supported function 9? on M. d , we have: 

lim / |y?^(k)| 2 dk — > 0, as R — > 00 
lim / |-^„| 2 (x)(ix — * 0, as R — > 00. 

^°>/|x|>R 

A practical sufficient condition is that ^(x, 0) is compactly supported and 
r/V'i/^x, 0) is square integrable with L 2 (M d )-norm bounded independently 
of 77. Then, we have the following convergence result [TBI [T5] : The Wigner 
transform W v (0, x, k) converges, after possible extraction of subsequences, 
in the space of distributions V'{R 2d ) to a Radon measure W (0,x, k), and 
moreover, we have 

/ t¥ (0,x,k)dxdk= lim / |-0,,| 2 (x, 0)dx. (6) 

In other words, the limiting Wigner transform captures all the energy of the 
incident wave field ip^ in the limit 77 — > 0. 

Kinetic Model. Upon using the Ito formula, we obtain that the average 
Wigner transform 

a v (t,x,k)=E{W v (t,x,k)}, (7) 
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solves the following kinetic equation 

da v f ~ dq 



k ■ V x a^ + Roa v = / R(k — q)a v (t, x, q) 



dt ' x " ' U ^ J Rd V ^^v^^^ (2?r)d , (g) 

a„(0,x,k) =W„(0, x,k), 

where we assume that ^(x, 0), whence a^(0,x, k), is deterministic; see e.g. 
[2] for the details of the derivation. We have defined R = R(0) and R(k) as 
the Fourier transform of i?(x), with the convention that 

R(k) = FR(k) = [ e- ikx i?(x)rfx. (9) 

Since -R(x) is a correlation function, R(k) is non-negative by Bochner's the- 
orem. For the rest of the paper, we assume that R(k) G PI L°°(]R' i ). 
Note that J R2d a v (t, x, k)rfxc?k is independent of time so that the total energy 
of the initial condition is preserved by the transport evolution. 

Scintillation. The validity of the kinetic model (jHJ) to describe the ensemble 
averaging of the phase space energy density of the wave field is trivial in the 
Ito-Schrodinger regime: the kinetic model (jHJ) is here exact for all 77 > 0, 
unlike what happens in other regimes of wave propagation [6], [TJ [19]. It 
remains however to understand how stable it is. In other words, how good 
an approximation is a v (t, x, k) of the random field W v (t, x, k). A natural 
object in the study of the statistical stability of W v is the following covariance 
function: 

J v (t, x, k, y, p) = E{W v (t, x, k) W v {t, y, p)} - E{W v (t, x, k)}E{W v {t, y, p)}. 

(10) 

We refer to this function as the scintillation function, in analogy to how stars 
are perceived to twinkle because the realization of the atmosphere changes 
in time. 

We shall see that the size of the scintillation function crucially depends 
on the smoothness of the initial conditions ^(x; 0) and a v (0, x, k) and on the 
support of the domain over which the energy density is averaged. The effect 
of the averaging will be quantified by measuring J v in appropriate (weak) 
norms. 

One of the main advantages of the Ito-Schrodinger regime of wave prop- 
agation is that J v (t, x, k, y, p) satisfies a closed form equation. Another ap- 
plication of the Ito formula [2] shows that J v is the solution of the following 
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kinetic equation: 

+ % + 2i?o - 0.2 ~ £ n ) J v = £ v a v ® Or,, (11) 
with vanishing initial conditions J v (0, x, k, y, p) = 0, where 

T 2 = k • V x + p ■ V y 
Q 2 J = J^^R(k-k')8(p-p') + R(p-p')8(k-k'))j(k',p')^^ 
*- r sr^ f A/ \ i (x - y),u , / , u u. du 

^ = 2^ / R ( U > e " e i e iK x > k + e ^>y>P+ e ilT)A7-w- 

. fi, 2 2 ( 27r ) 



(12) 

In the absence of the operator JC V , the variables (x, k) and (y, p) remain 
uncoupled in (fTTj) and the scintillation vanishes. Scintillation is created as 
the waves propagate through the random medium with a rate of creation 
proportional to K v a v <g> a,,. Notice that K v involves a highly oscillatory 
integral. Outside of the diagonal x = y, this oscillatory integral is small, 
whereas in the vicinity of the diagonal x = y, it is not. We thus observe that 
K v h is small when h is smooth and large when part of h is concentrated near 
x = y. 

Outline. The rest of the paper is structured as follows. The main results 
of the paper are summarized in section [2J We obtain estimates for J v in 
various norms, and in the specific case of initial conditions for a v of the form 
a n (0, x, k) = <5(x)/(k), show that rj~ l J rj converges to a measure J solving 
an explicit kinetic equation. Section [3] presents stability estimates for the 
scintillation operator K v defined in (1121) and for the kinetic equations (jHI) 
and (llip . The proof of the stability estimates for J v are given in section H] 
whereas the proof of convergence of rj J v when a^(0,x, k) = 5(x)/(k) is 
given in section [5) 



2 Main results 

Let ^(x, 0) be a sequence of rj— oscillatory, compact at infinity, functions 
uniformly bounded in L 2 (M. d ). This is the case of interest for us here, where 
we can define the Wigner transform (jSJ) and pass to the high frequency limit 
rj — > while still ensuring that energy is conserved as in §6§ . We are interested 
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in quantifying the statistical stability of the Wigner transform W n (t, x, k) and 
do so by analyzing the scintillation function J n defined in f llOp . 

We present two results. The first result proposes an upper bound for J v 
in different norms and for different initial conditions ^(x, 0). The second 
result analyzes the convergence properties of J v as rj — > for initial conditions 
of the form a^(0,x, k) = <5(x)/(k), which correspond to localized sources at 
position x = radiating energy smoothly in wavenumber k. In this context, 
we will show that J v is of order 0(rj) and will obtain the limit of rj~ l J v as 
Tj -> 0. 

Some typical initial conditions. Let us consider initial conditions ^{x., 0) 
oscillating at frequencies of order 7]~ l and with a spatial support of size r] a 
for < a < 1. The parameter a quantifies the macroscopic concentration of 
the initial condition. 

The simplest example is a modulated plane wave of the form: 

<'(x) = 4x(^)^, (13) 

where x(x) is a smooth compactly supported function on M. d . The direction of 
propagation is given by k . Note that the above sequence of initial conditions 
is indeed uniformly bounded in L 2 (R d ), compact at infinity, and ^-oscillatory. 
As another example of initial conditions, we consider 

where Jo is the zero-th order Bessel function of the first kind. Such an initial 
condition is supported in the Fourier domain in the vicinity of wavenumbers k 

(2) 

such that |k| = |k | so that ip^ emits radiation isotropically at wavenumber 
|k |; see [HI IS] for more details. We again verify that the above sequence of 
initial conditions is indeed uniformly bounded in L 2 (M. d ), compact at infinity, 

and ^-oscillatory. For this, we use that Jq{z) = cos(z — | ) + 0(z~ 3 ^ 2 ). 

Domain of measurements. For the above initial conditions for ip v , we are 

interested in the corresponding Wigner transform W v (t, x, k) and scintillation 
function J v . It turns out that J v is itself oscillatory so that its size depends on 
the scale at which it is measured. In order to capture this scale, we introduce 
a test function cp G 5(lR 2a! ), a fixed wavenumber ki G M d , and define 

1 / x k — ki \ 

^, sl , S2 (x,k) = -^- T - y ^^_,__j. (15) 
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We then denote by (•,•) the duality product iS'(IR n )-iS(IR n ) for n = 2d or 
n = Ad and want to quantify {W^, *p v , si ,s 2 )i the energy density averaged over 
a domain (in the phase space) of width rj si in space and 77 s2 in wavenumbers. 

By using the Chebyshev inequality, we obtain the following estimate on 
the probability that W v deviate from its ensemble average a v : 

P(lW),^, Sl , S2 )-KW)fe, S 2)l ><*) < -p{Jr,{t)iVr,,sx,S2®Vri,s x ,s 2 )- 

(16) 

Here, a <8> a(x, k, y, p) = a(x, k)a(y, p). In other words, when the above 
right-hand side converges to 0, then we find that (W v (t), (p v ,s 1 ,s 2 ) converges in 
probability to 0, which implies that W v (t) converges weakly and in probability 
to 0. The measured energy density is thus asymptotically statistically stable. 
A very relevant practical question pertains to the largest values of Si and S2 
that can be chosen so that the Wigner transform is still statistically stable in 
the limit 77 — ^ 0. We are now ready to state our main theorem on this issue. 

Bounds for the scintillation function. For any cp(x, k) G L 2 (M. 2d ), let 
jF x y9(u, k) and jF k </?(x, £) be the Fourier transforms of tp in the first variable 
only and in the second variable only, respectively. We also denote by a < b 
the inequality a < Cb, where C > is some universal constant. Then we 
have the following result: 

Theorem 2.1 Let ^(x, 0) be a sequence of functions uniformly bounded in 
L 2 (IR d ), compact at infinity, and rj- oscillatory. Let a^(0,x, k) be the cor- 
responding sequence of Wigner transforms given by ([5]). We assume that 
.T-xa^O) and !Fka v (0) are integrable functions and that 

11^x^(0, u,k)|| L i (K2d) < r]~ ad and H-F^O, x, £) ||li(rm) < V~ Pd , (17) 

for some a G M and (3 G R. Then we find that 

(Jr,(t),p v>Sl>S2 ® yw 2 > < r/-^ - 2 ^ ( v ^-*-^) A r, d ^ a -^). (18) 
Here, a Ab = min(a, b) and a V b = max(a, b). 

Of interest here is the following corollary: 

Corollary 2.2 Let ift v (0) be given by one of the expressions in ( |T3l) or ( 1141) . 
Then ( |T8l) holds with (3 — 1 — a.. 
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We can deduce the following results from the above corollary. In what fol- 
lows, we consider that averaging takes place over a large domain of wavenum- 
bers so that S2 = 0, as e.g., in spatial measurements of the physical energy 
density. 

Support of the sources. Let us assume that the spatial support of the 
domain of measurements is large so that s± = as well. Then we find that 



In other words, the scintillation is of order 0(r] d ) when a = 0, which cor- 
responds to a large support of the initial source term. This corresponds to 
the ideal case where the scintillation is smallest. In such a setting, we obtain 
that (W v — a v , ip) is of order 772 . This is the most stable situation. 

For a very narrow support of the initial source term comparable to the 
correlation length of the medium, namely when a = 1, we obtain that the 
scintillation is of order 0{rj) so that (W v — a v ,Lp) is now of order 772 . We thus 
obtain statistical stability of the energy density generated by a very localized 
source term whose radiation pattern in k is smooth, although the statistical 
instability is much larger than in the case a = 0. We know that for sources 
that are highly localized both in space and in wavenumbers, the scintillation 
does not converge to and the energy density is not asymptotically statisti- 
cally stable; see [2]. Such highly localized initial conditions would correspond 
to a choice a — [3 — 1 in Theorem 12.11 We will confirm in the next theorem 
that the order 0(rf) above is optimal. 

Small domain of measurements. Conversely, we can consider the case of 
a source term with a large support, which corresponds to a = 0, and a very 
small measurement domain. In this setting, we find that 



This means that the energy density becomes asymptotically statistically sta- 
ble as soon as it is measured over an area that is large compared to the 
correlation length of the medium. This is an optimal result of self-averaging 
as we cannot expect the energy density to be statistically stable point-wise, 
or when averaged over sub- wavelength domains. The above result, which is 
based on estimating K v in (lT2"j) in appropriate norms, improves on estimates 
obtained in 




(19) 



(20) 
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We can also consider intermediate situations where both the source and 
the measurement domain have small support. In that case, the optimal es- 
timate for the scintillation depends on whether a < S\ or si < a. These 
results are in fact optimal when the source term and the domain of measure- 
ments are located at the same place. Such a geometry explains why we do 
not obtain scintillation proportional to /^(i-a-si). We should obtain better 
estimates when the domain of measurements and the source term are not 
centered around the same point, though this cannot be inferred from our 
current results. 

Convergence of scintillation. Let us consider the case of initial conditions 
of the form (Tl3|) or (fT4|) with a = 1, i.e., for tightly localized source terms, in 
(transverse) dimension d > 2. The Wigner transform of such source terms 
converges in the limit 77 — * to a distribution of the form <5(x)/(k), where 
f(k) is a smooth function when x( x ) is smooth [T5]. We consider the kinetic 
equations with such initial conditions and obtain the following result. 

Theorem 2.3 Let J v be the solution of (fTTj) with the initial condition in (jSj) 
given by a^(0,x, k) = <5(x)/(k) for some smooth function /(k) in dimension 
d > 2. Then 7]~ l J ri {t) converges in the space of distributions uniformly in 
time to the limit J(t), which solves the following kinetic equation 

^ t +T 2 + 2R -Q 2 y=J°, (21) 

where 

f>(t) = e- 2Rot 5(x - tp)8(y - tq)2nx 

££(u)5(u • (p - q)) /(p - H) ( /(q _ |) _ /(q + |)) du. (22) 

The above theorem should be interpreted as follows. As time propagates, the 
transport ballistic part a°(t, x, k) = e~ i?0 *5(x — £k)/(k) creates some insta- 
bilities, which converge after appropriate scaling to J°(t). The scintillation 
thus generated is then transported by the transport equation (|2ip . We also 
observe that the error estimate of order 0{rj) in (Tl9l) with a = 1 is optimal. 



3 Functional setting and stability estimates 

In preparation for the proof of the theorems and the corollary presented in 
the preceding section, we prove here some stability results for the transport 
equations (jSJ) and ffTT]) and for the scintillation operator K n . 
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We denote by T the operator of Fourier transform with respect to all 
variables of the function on which it applies. For 1 < p < oo, we introduce 
X p as the subspace of tempered distributions in iS'(K 4,i ) such that 

\\hf Xp = sup f sup|^(u^,v,C)| p rfu<oo, (23) 
v,ceK d JR d £eR d 

for 1 < p < oo and 

IHk»= SU P |^/i(u,|,v,C)| < oo. (24) 



We also define Y p as the subspace of tempered distributions in <S'(IR 2d ) such 
that 

\\g\\ p Y = [ sup|^(u,OI P du<oc, (25) 

for 1 < p < oo and 

\\g\\ Yoo = sup \Th(u,£)\<oo. (26) 

Finally, we define Y as the subspace of tempered distributions in ^'(IR 2 ^) 
such that 

||flf||y = sup / |J 7 flf(u,^)|du< oo. (27) 



Morally (though this is inexact), the space X\ corresponds to scintillation 
functions that are integrable in one spatial variable (bounded in the corre- 
sponding dual variable v) and bounded in another spatial variable (integrable 
in the corresponding dual variable u). It is this boundedness that allows us 
to obtain the result fl20l) in the presence of small domains of measurements. 
In contrast, corresponds to scintillation functions that are integrable in 
both spatial variables (bounded in u and v) , which allows us to get the result 

(USD- 

The above spaces are well-adapted to the estimation of the scintillation 
operator K v . More precisely, we have the following result: 

Lemma 3.1 Assume that R £ L 1 (R d ) n L°°(R d ). Then for 1 < p < oo, 

(i) K v is bounded in X p and 

U^IUtXp) < 4||i?|| L i (IR d). (28) 
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(ii) Let /i G Y p and v G Y . Then 



\JC v fi®v\\ Xp < 4,r] d \\R\\ L ™Q SL d ) \\ii\\Y p \\v\\Y- 



(29) 



Proof. With obvious notation, we recast K v = e . e^e^K^ . Let h G X p . 
Then we have 

TK}i h = f e iw <^ + ^ ) R(w)J^h (u--,$,v + -,c) dw, 
so that using the Holder inequality with 1 



1 + 1 

p p 



hf Xp < sup / sup 



\R(w)Fh ( u 

V V 



w w 



^ II^IIlW) SU P / SU P / l^( W ) 

v,CeM d JR d £eM d ■ 



w w 

Th ( u ,£,vH ,C 

?7 77 



This proves (i). Let now /i := //<g)iA Upon performing the change of variables 
w — > rjw, we have 

TK^ii®v = rj d [ e ivw ^ ei * + ^ c) R{r]w)J 7 [i®v{u-w,€,v + w,C)dw, 



so that 
< rf sup / sup 



v. 



\R(r)(w — v))jF/x tg> v (v + u — w, £, w, £) |dw 



du, 



\l°°(r*)\\ u Wy su P / SU P / l-^Mv + u - w ^)| P |^( w .O |dwdu, 

v,CGK d </R d £eiR d ■ 

<^ll«ll^r»-JI*JI* 



< 77 d ||^| 



which concludes our proof. □ 

We need stability estimates for the kinetic equations. We start with the 
first kinetic equation: 

da 



dt 



+ p-V x a + R a = Qa + S, o(0, x, p) = a (x, p), (30) 
Qa(t,x,p) = (2n)- d [ R(p-p')a(t,x,p')dp', 
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with R Q := R(0), R G L 1 (R d ) n L°°(IR d ) and R non-negative. Then we have: 

Lemma 3.2 Assume that a$ £ Y p and 5 G L 1 ((0,T),Y p ) for some T > 
and 1 < p < oo. Then l[JU\) admits a unique solution in C°([0,T],Y p ) such 
that 

||a||c°([o,r],Yp) < \\o-o\\y p + ||5||li((o,t),y p )- (31) 

Let S = and let a°(t, x, p) := ao(x — £p,p)e _fioi 6e £/ie ballistic part of a. 
Then, assuming that T^a® £ L 1 (M 2d ) ; we /iave i/ie following estimate for all 
t > 0: 

\\(a-a°)(t,-)\\ Y <t^ d ! supl^ao^Ol^^^ll^kaollLX^). (32) 
Proof. The proof is a direct application of the integral formulation of 

a (t) = e- m g t a + [ e" Jlo ( t ->&_.g(a(s) + S(s))ds, 
Jo 

where Q t is the free transport semigroup given by 

(? t a(x,p) := a(x-tp,p). 

The operators Q and Q t are both continuous in Y p . Indeed, for ip G Y p , we 
have: 

so that 

< \\<p\W P , / 33 n 

\\Q<p\\y p < \\R\\l^^)\\v\\y p - 

Standard fixed point techniques then provide existence and uniqueness results 
for fl30|) . When 5 = 0, estimate (l3~Tj) follows from the maximum principle and 
the observation that ||a ||y p is a majorizing solution to fl30|) . When a = 0, 
( 13"Tj) is an application of the Gronwall lemma. 

For 5 = 0, we have the following Neumann series expansion in terms of 
multiple scattering: 

a n (t) = [ e- Roit - s) g t . s Qa n - 1 (s)ds, 
Jo 
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with the ballistic part a°(t, x, p) := e i?0 *ao(x — tp, p). By induction, we find 
the following expression for the Fourier transform of a n : 

ft PSl fS„-l 

JFa n (t,u,k) = e~ Rot / ■■■/ R(k+(t- Sl )u)--- 

Jo Jo Jo 

R(k + (s re _i — s n )u)jFa (u, k + tu)dsi ■ ■ ■ ds n . 

The change of variable k + tu — > £ yields 

e -Rot r 

\\a n (t,-)\\ Y < -jT^WRWl^) sup \Fa (v,$)\d£, 

g— Rot 

- n \ t d-n \\ R \\L^{M.<i)\\^ a o\\L^(R^). 

Summing over n > 1 gives the result. □ 

The last lemma deals with the fourth-order transport equation ( jTTl) : 

Lemma 3.3 Assume a £ X p and S £ L l ((0,T),X p ), for T > and 1 < 

p < oo. Then, the above system admits a unique solution in C°([0,T], X p ) 
such that: 

\\a\\c°({o,T],x p ) < ||«o||x p + IISIU^o.t),^)- (34) 
Proof. The result stems from the integral formulation of ( fTTI) given by 

a(i) = e- 2/?0 '^ 2 a + f e- 2R °^Q 2 t _ s [{Q 2 + JC v )a + S](s)ds, 
Jo 

where Qf is the semigroup defined as 

Qt a (^ Pi y, q) : = «( x - *p, p, y - *Qi q). 

From Lemma 13.11 we know K. v is continuous in X p , and so are Q 2 and Q 2 
since 

J"^ 2 V9 = !Fip(u, £ + tu,v, C + *v), 

= (i*(€) + i2(C))^(u ) €,v,C), 

for y> £ X p . Existence and uniqueness follow as before from standard fixed 
point theorems while estimate (]34p stems from separate applications of the 
maximum principle and the Gronwall lemma. □ 



14 



4 Estimates for the scintillation 

We are now ready to prove Theorem 12.11 and Corollary 12. 21 

Proof [Theorem I2.1j . According to Lemma I3.3[ the fourth-order trans- 
port equation (|2"T1) is stable in X p , so that we have the following estimate, 
uniformly on [0,T], 

W)lk< / \\)C v a v ®a v (s)\\ Xp ds. (35) 
Jo 

Provided that a v belongs to Y fl Y p , then K, v a v (g) a v is small in X p . Indeed, 
item (ii) of Lemma [3.11 yields for sG [0,T] and 1 < p < oo that: 

\\JC v a v ® a,, (s)||x p < 477^11^11^(^)11^(5)1^ ||a,,(s)||y p . 

First, we control the Y" norm by the Y\ norm since Y\ C Y. Lemma [3721 shows 
that the radiative transfer equation ([8]) is stable in Y r , for 1 < r < oo, so 
that we just need to estimate the initial condition a r; o(x, p) := 0^(0, x,p) in 
these Y r norms. Denoting by jF x a ??0 (u, p) the Fourier transform of a v o with 
respect to the spatial variable x only, we obtain, 



IKolk < / |^x«r,o( U )P)NP^ U ) 
IKolkoc < SU P / |^x«r,o(u,p)|rfp, 



so that the assumption of the theorem gives 

Kdk < c v - da . 

Moreover, defining ip v o(-) '■= ^»?(')0), we have the relation 

^ x a r ,o(u, p) = + | J ^VV^ (f " f ) ' 

from which it follows, using the Cauchy-Schwarz inequality, that 

\\a v o\\ Yoo < H-T^olli^) < C, 
where C is independent of rj. We have thus obtained that for all s G [0, T], 

H^a,®^)!^ < Crj d ^- a \ 
\\K Tl a v ®a ri {s)\\ Xl < C V d ^~ 2a \ 
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which yields by interpolation, for 1 < p < oo, 

||/C,a r/ ®a,( S )|U P <C ( 7 ? d ( 1 " (1+ ^ a ). (36) 

This induces a first estimate for J v , which is not optimal for initial conditions 
with small support when a — f3 > 0. The stability of the transport equation 
(IED in Y p is not sufficient to deal with such irregular initial conditions. Rather, 
we need to separate the ballistic part from the scattering part in the kinetic 
equation to obtain sharper estimates and thus introduce: 

a v (t, x, p) := a° v (t, x, p) + a°(t, x, p), 

where a°(£, x, p) = e _i?0 'a^o(x — tp, p) is the ballistic part and a s v satisfies 

da s 

-Qf + P- V x a; + Ro a s v = Qa s n + Qa% a s v {t = 0, x, p) = 0. 

Since the Fourier transform of a° is given by e- Roi J r a ?? o(u,k+tu), its Y norm 
can be estimated for t £ (0, T] as: 



\<(t)\\v < 



sup / iJFa^u, k + tu)\ du < -g sup / |jFa r? o(t~ (k — z),z)|dz 

If c 

< — d / |^kSo(x,z)| rfzdx < — d vT dp 

t J R 2d I 



Now, Lemma [3.21 and estimate fl32l) imply that: 

so that the time singularity of a s v is weaker than that of a°. Thus, for 1 < 

p < oo, 



Ka v <g> a v (s)\\ Xp < ^M~H (s- d + s l ~ d ) < 



s~ d . 



For short times, we then use estimate (1361) since it is independent of s and 
for longer times, we use the above estimate. We thus write: 

\\)C n a n ®a v (s)\\ Xp = I(s < ^o(^)) H^S^S ( s ) \\x P + ^ > t (rj)) \\K v a v ®a v {s)\\ Xp 
so that, for t £ [0,T], we have 
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Setting t Q (rj) = r] a 13 when a > [3 above and using (155]) and (1551) . we find, for 
t E [0,T], that 



I J, 



< (7r ? a! ( 1 ~( 1+ |) Q )+( Q -/ 3 ) v0 . 



We conclude by using the Parseval-Plancherel equality which yields, for t G 

[o,n 



< (2tt) d || JqHx, 1 1 fywi, s 2 1 1 L 1 



I^Vu,«i,«(' J P)l d P 



with 1 = i + 4. It remains to verify the scaling properties: 



||-^<A7,si,s 2 |Ui(]R2<i) - ^(5!+^) ll-^V 9 l|Li(K 2£i ) 
l^r,, S1 , S2 (-,P)|^P 



d( Sl (l-I) +S2 ) 



|^(-,p)|dp 



We conclude the proof of the theorem by choosing p = oo or p = 1 in the 
above estimates. □ 

Proof [Corollary I2.2j . We simply need to estimate T^a^ and T^a^ in 
L 1 (M d ). Since 

^*So(u, P) = + ^ ^ (| - ^) ' 

J'kSoCx, k) = VryO (x + |kj (x - ~kj , 



it follows that: 



( 5 + H j T ^ _ H 



/ |^x«r,o( U )P)l rfurf P = / 

= / |^r?0 ( U ) Fl/)lfl (p)Uurfp = ||^r,o||ii(Rd) < 
J R 2d 

l^oC^pjidxdp^^-i^ollii^) < cv d{1 ~ 0) - 



dudp, 



da 



It suffices to set (3 — 1 — a in Theorem 12.11 to conclude the proof of the 
corollary. □ 
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5 Convergence of the scintillation 

We now prove the announced convergence result. We first observe that the 
existence results obtained in Lemmas 13 . 21 and 1 3 . 31 hold when the spaces Y p and 
X p are replaced by the spaces of bounded measures A4(R 2d ) and J\A(M. M ), 
respectively or by the spaces of continuous functions C°(JBL 2d ) and C°(M. M ), 
respectively. We recall that d > 2 here. 

Proof [Theorem 12.3] . The scintillation function satisfies the following 
transport equation in integral form 

J v (t) = [ e- 2R * it - s) gl s {Q 2 + K v )J v {s)ds+ f e- 2R ^- s) g 2 t _ s K v a®a(s)ds. 
Jo Jo 

(37) 

We recast this, with obvious notation, as 

oo 

T _ rp j , tO t_ \ ^ rpk jO 

V It] 7] ' rj ' °v ~ / 4 ± 2n°n- 

k=0 

We denote by T 2 the formal limit operator of T 2j? defined as 

T t f(t)= f e- 2R ^g 2 _ s Q 2 f(s)ds. (38) 
Jo 

The source contribution. We verify that 



K v a ® a(s) - K v a° ® a°(s)|| Xoo <-^Al, (39) 



where the ballistic part is given by 



a°(t, x, k) = e"^(x - *k)/(k) = e-^5(k - 

Indeed, we know from Lemma 13.11 that 

\\K v (a - a ) ® a(s)\\ Xoo < r, d \\a - a°||y||a|| yoo , 
and from (1321) in Lemma [3.21 that 



la — a°||y 



< t 1 " 1 [ sup |^ao(v,£)|d€ < t 1 ^ / 
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That Ijally^ is bounded comes from the stability of the transport equation in 
established in Lemma T3.21 The term K n a° (8) (a — a ) is treated similarly. 
Let us define 



(40) 



joo(t) = / e- 2Ro ^gl s K v a ®a°(s)ds, 
t 

= e~ 2Rot [ gl s K v g 2 s (a ® a )ds. 
Jo 

We find that 

\\J v (t)-J v (t)\\ Xoo < V ^^ V . (41) 
Indeed, we deduce from (1391) and the stability of Qf in C^X^) that 

- J?(t)\\ Xoa < l\4=l A < *o + V d tl~ d < »A 
Jo s 

d 

for to = V 1 ' 1 • Up to a smaller-order error term in the space of distributions, 
we may thus replace J° by J°° in the sequel since the transport equation (1371) 
is stable in X^. Now, calculations with K v replaced by K~ l ~ l show that 

e 2 ^J°° n (t) = f j RH e^ [{x - ( *- s)p) - (y ~ ( *" s)q)] (5(x-tp + S f) 

*(y-*q+sf)/(p-f)/(q-f)duds 
= J J R(u) e^- (p_q) 5(x-tp + sf) 

*(y-*q+«f)/(p-f)/(q-f)*«fa, 

_^ 

e 2flot J° n (t) = r] J n J R(u) e isu ( p - q )<5(x-tp + ?7sf) 

5(y-tq + t7flf)/(p-H)/(q-f)du£fa. 
Upon sending 77 — > 0, we find in the limit that 

7 00 (t) 

r/^0 7] 

e -MoM(x - tp)<5(y - tq)vr / i?(u)5(u ■ (p - q))/(p - H)/(q _ ") du 

+ze- 2 ^^(x - tp)*(y - tq) p.v. / i?(u) — J -/(p - H)/( q - H) du , 

jRrf u • (p - qj ^ ^ 
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in the space of bounded measures A4(R 4d ). After accounting for all four 
terms in the definition of K v and using the fact that R(u) = R(—u), we find 
that the limit of ?7 -1 </°(£) is given by: 

J°(£) = e - 2Ho *5(x - tp)6{y - tq) x 

2tt J R(u)5{u • (p - q)) /(p - H) ( /(q _ «) _ /(q + H)) du. 

This gives us the source term (122!) in the transport equation (j2Tj) . 

Kinetic equation for the scintillation. We have shown that ri" 1 ,]® con- 
verged to J°. It remains to obtain convergence of the whole sequence r]~ l J v . 
Let x, p, y, q) be a a smooth function on [0, T] x IR 4d . Then we have by 
integration on the latter space that 



(J,,0) = (T 2 ,J^) + (J° 0), (42) 



and equivalently that 



(J,,0) = (J,,T^) + (J° 0), (43) 



with 



T 2 >( S ) = / T e- 2i ^(g; + K;)Q£ s <i>{t)dt. (44) 



We have shown that the difference between the source terms r\ ' J2 and 
J° converges to as a distribution and has a negligible effect on r/" 1 ^. So 
we can replace the initial condition for the error term by J° and look at the 
problem 

where J v is now of order 0(1). We observe that 

J°(t) = e - 2Rot 5(x - tp)5(y - tq)H(p, q), 

where H(p, q) is a smooth function. 

Let now J* = J v — J° be the solution of 



We recall that 



T 2n J(t) = f e- m ^gl s {Q2 + K 11 )J{s)ds : 
Jo 
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so that 

T2rjJ° = T 2 J° + J^, 

where J 2 is given by a bounded operator in .M(IR 4d ) applied to K V J°. The 
latter is given by 

<5(x - tp)S(y - tq)e~ 2 ^ [ R(vL)e iE ^- u H(p - H) #( q - ^)du, 

plus similar contributions. Because H is a smooth function, this term con- 
verges to in Jvt(M. 4d ) as rj — > 0. This shows that J 2 converges to as 
rj -»• 0. 

The other contribution, T 2 J°, involves a bounded operator applied to 
Q 2 J°, which is equal to 

Q 2 J°(i)(x, P , y, q) = e- 2 ^i?(p - *)i?(q - f , f )• (45) 

For /, whence if, and i? sufficiently smooth, the above function is bounded 
in C°(IR 4d ). The function is not bounded uniformly in time, however, and we 
split the contribution J°(t) into J$(t) = J°X(o,s)(t) an d J°X(s,T)(t), which we 
still denote by J°(t). The source term T 2 J® generates a small contribution, 
which goes to as 5 goes to in the sense of distributions since the term 
in (f45j) is bounded in e.g. L 1 (R 4d ) uniformly in time so that after time 
integration in (I38j) . the contribution is bounded by 0(5) — > 0. The remaining 
contribution is bounded in the uniform norm uniformly in time with bound 
inversely proportional to 5 2d . 

We now have a problem of the form 

J\ — T 2v J^ + T 2 J°, 

where T 2 J° is uniformly bounded in the uniform norm by 0(5~ 2d ). Weakly, 
this means that 

(j, 1 ,0) = (j, 1 ,t 2 ;0) + (t 2 j°,0), 

where (f)(t, x, p, y, q) is a smooth function. The solution J3 is bounded in 
C°(M 4d ) uniformly in 7] by stability of the fourth-order transport equation in 
the uniform norm. There is therefore a subsequence that converges weak * 
in L°°(R 4d ) to a limit J 1 e L°°(M 4d ). 
Let us decompose T 2ri as: 

T 2v J(t) = f e- 2R °^gl s (Q 2 + K v )J(s)ds = T 2 J(t) + S 2r , J(t), 
Jo 



21 



where 

s 2v J(t) = f e-^-^gl^j^ds. 

Jo 

We choose <j) sufficiently smooth so that S% (f) goes to strongly in L 1 (M. M ). 
As a consequence, (J*, S^cj)) goes to with rj so that, in the limit, we have 

(j\0) = (j\T^) + (T 2 J°,0). 

The above convergence to the limiting transport equation holds for every cut- 
off 5. Thus, by stability of the limiting transport equation, we can remove 
the cut-off in 5 and obtain that 

J 1 = T 2 J 1 +T 2 J°, 

weakly in the space of distributions. The above integral equation admits a 
unique solution, which shows that the whole sequence rj converges to J 
solution of: 

J = T 2 J+J°. 

This completes the proof of Theorem 12.31 □ 
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